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ABSTRACT 

We exploit a previous computation of the self-mass-squared from quantum 
gravity to include quantum corrections to the scalar evolution equation. The 
plane wave mode functions are shown to receive no significant one loop cor- 
rections at late times. This result probably applies as well to the inflaton 
of scalar-driven inflation. If so, there is no significant correction to the i-ptp 
correlator that plays a crucial role in computations of the power spectrum. 
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1 Introduction 



Massless, minimally coupled scalars and gravitons are unique in possessing 
zero mass without classical conformal invariance. Simple arguments based 
on the energy-time uncertainty principle suggest that quantum effects from 
these particles should be vastly enhanced during primordial inflation [1]. 
The lowest order effect is that inflation rips long wavelength virtual quanta 
of these types out of the vacuum. That is the origin of the observed scalar 
perturbations ^ and of the potentially observable tensor perturbations [3]. 

Higher order effects derive from the interactions of the sea of long wave- 
length virtual scalars and gravitons, both with themselves and with other 
particles. Loop corrections to the primordial power spectra are enhanced by 
logarithms of the ratio of the scale factor to its value at first horizon crossing 
[ll[5]. For any mode whose spatial variation we can resolve today these sec- 
ular enhancements cannot overcome the minuscule loop counting parameter 
of GH"^ < 10^^^. Exceptions have been suggested [SI 13 IB]- However, much 
larger loop corrections (to other things) can be obtained either from interac- 
tions with larger loop counting parameters, or by studying things we would 
perceive as spatially constant such as the vacuum energy or particle masses. 

Explicit computations have been made on de Sitter background in five 
different models which involve either scalars or gravitons: 

1. For a massless, minimally coupled scalar with a quartic self- interaction 
— the two loop expectation value of the stress tensor [H] and the two 
loop corrections to the scalar mode function [TUl ITT] . 

2. For a massless, minimally coupled scalar which is Yukawa-coupled to 
a massless fermion — the one loop corrections to the fermion [121 [13] 
and scalar [T3] mode functions, and the two loop vertex function [T5] . 

3. For massless, minimally coupled scalar quantum electrodynamics — 
the one loop correction to the photon [161 [H] a^iid scalar [T8| mode 
functions, the two loop corrections to the coincident scalar [19] and 
field strength [2^ bilinears and to the stress tensor [20j. 

4. For general relativity — the one loop self-energy [21] and the expecta- 
tion value of the graviton field at one [22l [23t [2l] and two loops [25] . 

5. For general relativity coupled to a massless fermion — the one loop 
corrected fermion mode functions [25] . 
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The leading secular effects of the first three models agree precisely with the 
stochastic formalism of Starobinskii [271I2S1I2S], and the series of these leading 
effects at all orders can be resummed [30l [151 EI] to give nonperturbative 
predictions. No stochastic results are available for the last two models. 

It is natural to extend these studies by combining general relativity with 
a massless, minimally coupled scalar to probe the effects of gravitons on 
scalars. Such an investigation has great phenomenological interest because 
the inflaton potential of scalar-driven inflation is so flat that inflaton mode 
functions are effectively those of a massless, minimally coupled scalar. In- 
deed, it is standard to compute the power spectrum of scalar perturbations 
by setting the scalar-scalar correlator to its value for a massless, minimally 
coupled scalar on de Sitter background [321 E3]- The subject of quantum 
gravitational corrections to the scalar effective potential has a long history 
[3H |35] but we will here look at all one loop corrections to the linearized, 
effective scalar field equation, including corrections to the derivative terms 
and also the fully nonlocal corrections. This is what one must do to fix the 
scalar field strength in addition to its mass. It is worth noting that one loop 
quantum gravitational corrections induce a time dependent field strength 
renormalization for massless fermions which eventually becomes nonpertur- 
batively large |26j . 

Our result is that the scalar mode functions experience no significant 
corrections at one loop. We prove this by solving the linearized, effective 
scalar field equation. Section 2 gives a brief review of our previous one 
loop computation of the scalar self-mass- squared ^36j, which represents the 
quantum correction to the linearized effective field equation. In section 3 
we actually solve the equation in the relevant regime of late times. Our 
conclusions comprise section 4. 

2 Our Previous Calculation 

The linearized effective scalar field equation consists of the classical term 
minus the scalar self-mass-squared M^{x;x') integrated against the scalar, 

^^(^^/^g'"'^Mx)) - I d*xM\x-x')^{x') = . (1) 
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We computed M^(x; x') at one loop order using the bare Lagrangian of Ein- 
stein + Scalar, 

C ^ -\d,^ d.^g^^^g + ^ {r - {D-2)a) . (2) 

Here A = {D — 1)H'^ is the cosmological constant and G is Newton's constant. 
We worked on the open conformal coordinate patch of D-dimensional de 
Sitter space, with the graviton field hf^i, defined as, 

giMi^iv, ^) = o^iv) (vpiv + f^h^^{r], where a{r]) = --^ . (3) 

Here 77^,^ is the Minkowski metric and = IGvrG is the loop counting param- 
eter of quantum gravity. The graviton propagator was computed by adding 
a non-de Sitter invariant gauge fixing term [STJ [38] , 

Cgf = -^V^g'^'F.F,, (4) 

= v'''{h^P,a-lh,^,^ + {D-2)aHh^X) ■ (5) 

The computation of M^(x;x') was made using dimensional regulariza- 
tion. We subtracted off the one loop divergences using three dimension six 
counterterms, 

AC = ^a,K^^,^,ip,p^g^''gP''^ 

a2K^R — a^K^R ^. — 

y/-g - ^j^ij^ V,t^,j9W-9 ■ (6) 



2D{D-l) ' ^ 2D{D-l) 

In these expressions a comma denotes ordinary differentiation while a semi- 
colon stands for the covariant derivative. The noninvariant counterterm pro- 
portional to as was made necessary by our non-de Sitter invariant gauge 
fixing function ([5]). Because gravity is not perturbatively renormalizable, 
only the divergent parts of the coefficients are fixed, 

ai = + Acti , (7) 
"2 = — TTi^ +^"2, (8) 



(47r)T t 3 

H'^-^ r 4 58 , , 





Coefficient of 


/ /\4i — \?, 1 tt2 

[aa j / H 


In a; 
3x 


1 l\A\ — \2 

[aa] VA 


26 Inx 1 38 1 2 i o 1 

" + 1^ — 6 m X — 18 m X 

3x 3x 


/ /\4 — 1 

\aa) H VA 


6 In X 1 4 /II 

X X 


/ /\4 rr4 

(aa ) ii 


4 In X 1 18 -1 /~io 1 1 o/~i 

^-^^ + — — 108 mx — 120 

X X 


/ /\3/ 2 1 /2M — ^3 1 tt2 

[aayia^ + a jD^/i/ 


In X 
6x 


/ /\3/' 2 1 /2M — \2 

[aayia + a )□ 


In X 1 1 
3x 6x 


(aa J (a + a j/i U 


21nx 1 5 101^,^ 

^ h - — 18 mx 

iX X 


(aa')^(a2 + a'^)H^ 


4 In X 32 ^4 
3x 3x 


{aa'fH^W^ 


2inx 16+841^3; 48xlnx + 96x 

3x X 




In X 
3x 


(aa')'(a' + a'^)H^V^ 


^''"^ + ^^^ 121nx + 48xlnx + 12x 

3x 2x 


{aa'fia^ + a'^)V^n 


i7inx 49^41^2 ^Qj^ -^2a; In X 

6x 6x 


(aa')'V^ 


y In X — 24x In X + 24x^ In x — 36x^ 



Table 1: Nonlocal Contributions to — ■jM^(x;x'). Here x = | and y{x;x') is 
defined in equation (IT^ . 



The finite parts Aoj are arbitrary. 

After renormalization and taking the unregulated limit {D = 4), our 
result is [SS] . 

-iM2(x;x') = m2a^(Aaitf+A«2n+Aa3^)5^(x-x')+ (Table [ij) . (10) 

The Laplacian and the d'Alembertian □ are, 

V^ = dA and n = \d''(a%) =-^d',--Hdo + \v^ . (11) 

a* V / a 

We should also define the de Sitter invariant length function that appears in 
Table [B 

y{x] x') = a{r])a{ri')^^\\x — x'\\'^ — {\r]—r]'\—i5)'^^ . (12) 
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3 Effective Mode Equation 



This is the heart of the paper. We begin by clarifying what is meant by the 
effective mode equation, then we explain the restricted sense in which we 
solve it. Finally we work out the contributions from the local counterterms 
in ( JTOl) and from the nonlocal terms of Table [H 

3.1 The Schwinger-Keldysh Formalism 

We seek to find plane wave mode solutions to "the effective field equations." 
Although the quoted phrase is common parlance, it is nonetheless ambiguous 
because there are different sorts of effective field equations whose solutions 
mean different things in terms of the unique canonical operator formalism. 
Introductory courses in quantum field theory typically concern the in-out 
effective field equations. The plane wave mode solutions of these equations 
give in-out matrix elements of commutators of the full field with a tree order 
creation operator of the in vacuum. 



This quantity is of great interest for fiat space scattering problems but it has 
little relevance to cosmology where there may be an initial singularity and 
where particle production precludes the in vacuum from evolving to the out 
vacuum. 

The more interesting cosmological experiment is to release the universe 
from a prepared state at finite time and let it evolve as it will. The mode 
solutions of interest to this experiment are the expectation values of com- 
mutators of the full field with the tree order creation operator of the initial 
vacuum. 



Keldysh formalism f2E[ [TB] . This is a covariant covariant extension of Feyn- 
man diagrams which produces true expectation values instead of in-out ma- 
trix elements [391 SOl El US]- Because there are excellent reviews on this 
subject [m HH HSl HH], we will confine ourselves to explaining how to use 
the formalism. 

The chief difference between the Schwinger-Keldysh and in-out formal- 
isms is that the endpoints of particle lines have a ± polarity. Therefore, every 




(13) 
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propagator iA{x\x') of the in-out formalism gives rise to four Schwinger- 
Keldysh propagators: i/S.^ 



and iA, 



Each of these propagators can be obtained by making simple changes to the 
Feynman propagator. For our model, the Feynman propagators of the scalar 
and graviton happen to depend upon the length function y{x; x') defined in 
expression f[T^ . and also upon the two scale factors. The four polarities 
derive from making the following substitutions for y{x; x'): 



i/S.+^{x] x') 



zA+_(a;; x') 



iA^+{x;x') : 



iA (x:x') : 



= a{ri)a{ri') — — (jr]— r^'l — 2(5)^ , (15) 
= a{vi)a{rj') — — {rj — rj'+i5Y , (16) 



y y++{x;x') 
y y+-ix;x') 

y y^+{x;x') = a{r])a{r]') \\x — x'\\'^ — {r] — r]' — i6Y , (17) 

y^y {x; x') = a{r])a{r]') \\x — x\\'^ — {\'r] — 'r]'\+iSY . (18) 



Vertices in the Schwinger-Keldysh formalism either have all + lines or all — 
lines. The + vertex is identical to that of the in-out formalism, whereas the 
— vertex is its conjugate. 

Because any external line can be either + or — , each N-point one particle 
irreducible (IPI) function of the in-out formalism gives rise to 2^ IPI func- 
tions in the Schwinger-Keldysh formalism. The Schwinger-Keldysh effective 
action is the generating functional of these IPI functions. We can express it 
in terms of fields <y9+, to access the -|- lines, and to access the — lines. 



d'^x / d'^x' 



+0(^i), (19) 



ip+{x)Ml^{x; x')(f+{x') + if+{x)Ml_{x;x')(f^{x') 
+ip^{x)M^^{x; x')ip+{x') + ip-{x)M^^{x; x')ip-{x') 

where S[(p] is the classical scalar action. 

At the order we are working, —iM^^{x; x') is the same as the in-out self- 
mass-squared. We can therefore read it off from (ITO!) . We get —iM^_{x; x') 
by dropping the delta function terms, reversing the sign and replacing y{x; x') 
by ?/+_(x;x') in Tabled! The other two IPI 2-point functions derive from 
conjugating these two. 



[x\ X 



-iMl^{x; x' 



-%Mt^{x\ x' 



(x: X 



(20) 
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To get the Schwinger-Keldysh effective field equations one varies the ac- 
tion with respect to the field of either polarity, then sets the two polarities 
equal to At linearized level this gives. 



Here rji = —1/H is the initial (conformal) time at which the universe is 
released in free Bunch-Davies vacuum. One can see from relations (fT5|) and 
(fTB]l . and from the extra conjugated vertex in M+_(x; x'), that the bracketed 
term vanishes for t]' > rj. The fact that y+_{x\x') is the complex conjugate 
of x') for rj' < rj means that the bracketed term is real. One also sees 

that it must involve the imaginary part of at least one propagator, which 
means the only net effect comes from points x''^ on or inside the past light- 
cone of x'^. Hence the Schwinger-Keldysh effective field equations are real 
and causal, unlike those of the in-out formalism. 

3.2 Restrictions on Our Solution 

Two limitations on our knowledge impose important restrictions on the sense 
in which we can solve (1211) : 

1. We only know the scalar self- mass-squared at one loop order; and 

2. We took the initial state to be free, Bunch-Davies vacuum. 

The first limitation means we must solve (!2T|) perturbatively. The full scalar 
self-mass-squared can be expanded in powers of the loop-counting parameter 



a 




(21) 



= levrG, 



oo 




(22) 



A similar expansion applies for plane wave solutions to fl2Tl) . 



oo 




(23) 



To make $(x; fc) agree with (|T4|) we must normalize the tree order solution 
appropriately, 

$o(r/, k) = uiri, k) = (l - ^) exp [^1 . (24) 
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The i > 1 solutions obey, 



a' 



+ 2Hdo + P ^e{r], k) = f dr]' fd^x'Ml{x; x')^e-k{v', A;)e*(^"-^1 



(25) 

We know only Ai1{x; x') so the sole correction we can compute is $i(?7, k). 

The second limitation means it only makes sense to solve for $i(?7, k) at 
late times, i.e., as — >■ 0~. Interactions result in important corrections to free 
vacuum on a flat background and it is unthinkable that this does not happen 
as well for de Sitter background. In the Schwinger-Keldysh formalism these 
corrections would correspond to vertices on the initial value surface |22l [18] . 
In the in-out formalism the free vacuum is automatically corrected by time 
evolution. One can follow the progress of this in the Schwinger-Keldysh 
formalism by isolating terms that decay with increasing time after the release 
of the initial state. For example, the two loop expectation value of the stress 
tensor of a massless, minimally coupled scalar with a quartic self-interaction 
gives the following energy density and pressure [9], 

Aif^ f 1 2/ X 1 1 >^ (n + 2)a-"-n , , 

(27r)4\ 8 ^ W2 ^ ^ 24^ (n + l)2 J v ; v ; 

We suspect that the (separately conserved) terms which fall like powers of 
1/a can be absorbed into an order A correction of the initial state. On the 
other hand, the terms which grow like powers of ln(a) represent the effect of 
inflationary particle production (in this case, of scalars) pushing the field up 
its quartic potential. 

Because we have not worked out the order k and corrections to Bunch- 
Davies vacuum, we cannot trust corrections to $i(?7. A;) that fall off at late 
times relative to the tree order solution uiji^ k), 



P 




«(^' ^) = ^ 1 + 7777^ + ^TTT:?!^ + O • (28) 



To understand what this means, it is best to convert equation ( ]25ll for $1 
from conformal time rj to comoving time t = — ln{—Hi])/H, 



$1 = rdr]'u{r]', k) jd^x'Ml{x- x')e'^-(^'-^\ (29) 
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At late times we can drop the factor of k'^ /a? on the left hand side. Given a 
putative form for the late time behavior of the right hand side it is easy to 
infer the leading late time behavior of $i, for example, 

r.h.s— .In(a) =^ _A_Z , (30) 

r.h.s^l ^ ^1--^, (31) 

r.h.s-.i^ ^ Constant - , (32) 

a 2H a 

r.h.s — > - =^ $1 — > Constant -— . (33) 

Because any constant contribution to $i could be absorbed into a field 
strength renormalization, the only effects which can be distinguished from 
corrections to the initial state derive from contributions to the right hand 
side that fall off no faster than l/ln(a). It turns out that inverse powers of 
ln(a) cannot occur, so the practical dividing line is between contributions to 
the right hand side which grow or approach a nonzero constant and those 
which fall off. 



3.3 Local Corrections 

Because the scalar d'Alembertian annihilates the tree order solution, only 
the third, noncovariant counterterm makes any contribution to (!29!l . 



J \ J 



XM(V,A;)e*'=-(^'-^^) = -Aag— M(r/,fc) . (34) 



This term rapidly redshifts to zero and we see from the preceding discussion 
that it cannot make a significant contribution to $i at late times. 

3.4 Nonlocal Corrections 

Although Table [1] might seem to present a bewildering variety of nonlocal 
contributions to (!29l) . a series of seven straightforward steps suffices to eval- 
uate each one: 
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^ °Vf^U+3, (35) 



1. Eliminate any factors oil/y using the identities, 

- = -M- 

2. Extract the factors of □ and from the integration over x''^ using the 
identities, 



□ 



2 . d fc2 
a TT^ + 4a— + 



da"^ da a?H'^- 



(37) 



V' -k' . (38) 

3. Combine the ++ and H — terms to extract a factor of i and make 
causahty manifest, 

In(^) -ln(^) =2vrz^(Ar/-Ax) , (39) 
ln2(^^) _ln2(^^)=47r2^(^A77-Ax)ln(^aa'i7'(Ar/'-Ax2)) . (40) 

Here we define a = 0(77), a' = ci{ri'), Arj = rj — rj' and Ax = ||x — a;'||. 
Note that any positive powers of y become, 

y+i — > -^aaH'^{Ar]'^-Ax'^) . (41) 

4. Make the change of variables f = x' — x, perform the angular integra- 
tions and make the further change of variable r = Arj ■ z, 

J £x'e(^Ari-Ax^F(^^aa'H\Ari'^-Ax^)y'''<'^'-'^'^ 

= 47re(Ar/) p\r r^pi^^aa' H\Ar^^ -r^))^^^^^ , (42) 

--(-'')-'''r-^^<-'(^j)'^(^))^-(«) 

5. Reduce the z integration to a combination of elementary functions and 
sine and cosine integrals [H]. 
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6. Make the change of variables a' = —l/Hrj', expand the integrand and 
perform the integration over a'. 



7. Act any derivatives with respect to a. 

Much of the labor involved in implementing these steps derives from the 
spacetime dependence of the zeroth order solution, k)e'^^'^' . For example, 
one can see from fH51) that only elementary functions would result from the z 
integration if the zeroth order solution were constant. In fact it is completely 
justified to make this simplification. 



1 





dr]'ld^x'Ml{x- x')u{r]', k)e"'-^'^'-^^ 

u{0,k) 



a4 





dri'ld^x'Mlix;x'); (44) 



Vi 



To see why, note from expressions fl28l) and fHSj) that the deviation of the 
zeroth order solution from u{0, k) introduces at least two factors of 1/a' or 
1/a. Because the constant mode function can at best result in powers of 
ln(a), and because it makes no physical sense to retain contributions to (129!) 
which fall off, we can immediately discard terms which acquire an extra factor 
of 1/a. Extra factors of 1/a' effectively restrict the integration to early times, 
which again causes the net result to fall off at late times. 

One consequence of the simplification (jUj) is that we can neglect any 
contribution from Table [1] which contains a factor of V^. We therefore need 
only compute terms of the form, 

-''("•*)^(^)7*''«-/.v{/(l^^)-/(^)} , ,4, 

where the constant K takes the values of 3, 4 and 5, and the functions f{x) 
are 1/x, \nx/x, lux and In^x. The action of the d'Alembertian derives from 
setting /c = in (1371) . The integrations for i^' = 3 (given in Table [2]) and 
K = 4 (given in Table [3]) were worked out in a previous paper [19]. Applying 
the same technique — which is just the 7-step procedure given above for 
k = — gives the results for K = 5 in Table HI 

It remains just to act the d'Alembertians and sum the results for each of 
the contributions from Table [H That is done in Table [51 Although individual 
contributions can grow as fast as In^(a), all growing or even finite terms 
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1 




\n{x) 

X 


ln2(a) 1 21n(a) 3 , tt^ , ^fln{a)\ 
2a ' a a ' 3a ' a2 ^ 


ln{x) 


l-'-^ + h + oi^) 


ln2(x) 






Table 2: Integrals with a'^. 




fir) 


ill* r 71 ' '1 r ^ ^/v^-m 
10^2 X j (1 .V a /(-'^ )} 




1 

X 






\n{x) 

X 






\n{x) 


lln(a)-§ + 0(i) 




ln^{x) 


ilnV)-iln(a) + |-^ + 0(i^^) 






Table 3: Integrals with a'^. 






-igx/dVan/m-/(V)} 




X 






ln(x) 

X 






\n{x) 


^a- | + 0(i) 




ln2(j;) 


13 1 , / X 1 5 1 /n (' ^'^('^) "1 
144« 3ln(a)+9 + U( ^ j 





Table 4: Integrals with a'^. 
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cancel. It follows that the right hand side of fl29p falls off at least as fast as 
1/a times powers of ln(a). Hence there are no significant corrections to the 
mode function at one loop order. 

4 Discussion 

We have solved the one loop-corrected, linearized effective field equation 
for a massless, minimally coupled scalar interacting with gravity during de 
Sitter inflation. Unlike previous analyses of the scalar effective potential [35] . 
our technique uses the full self-mass-squared, including corrections to the 
derivative terms and nonlocal corrections. It should therefore be sensitive 
not only to the scalar's mass but also to its field strength. We find no 
significant corrections of any sort at one loop order. 

Physically our result means that the sea of infrared gravitons produced by 
inflation has little effect on the scalar. Although our scalar is a spectator to A- 
driven inflation, typical scalar inflaton potentials are so flat that the inflaton 
is also unlikely to suffer significant one loop corrections from quantum gravity. 
If so, it is consistent to use the tree order scalar mode functions to compute 
the power spectrum of cosmological perturbations. 

There are good reasons why the scalar should not acquire a mass from 
quantum gravity |i35jj but there seems to be no reason why its field strength 
cannot suffer a time dependent renormalization of the type experienced by 
massless fermions [26]. Weinberg's result for the power spectrum allows for 
ln(a) corrections [1], and they do arise from individual terms in Table [S] 
However, there is no net correction of this form at one loop. 

In the end, our null result may not be as surprising as it seems. It is very 
simple to show that the scalar self-mass-squared vanishes at one loop order on 
a flat background [36] . It is also wrong to think of the scalar as weaker than 
the fermion; what is weaker is the fractional one loop correction. Although 
massless fermions experience a time-dependent field strength renormaliza- 
tion, this represents a ln(a) enhancement of tree order mode functions which 

3 

fall like l/aa. By contrast, the tree order scalar mode functions approach a 
nonzero constant. Indeed, it is the fact that derivatives of this constant van- 
ish which makes the fractional correction of the nhdipdip coupling so small. 
The sea of infrared gravitons is present, but it can only couple, at this order, 
to the scalar's stress-energy, and the stress-energy of a single scalar redshifts 
to zero at late times. 
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External 0])era.tor x f(x) 


Coefficient of (/(O./r) x ^f", 


{aa'ynyn^ X 










(aa')^n^ X f 





(aa')^n^ X -6\n^x 


-18 


(aa')^n^ X -181nx 





(aaO^i/^n X 6'"^=^ 

^ ' a; 





(aa')^i/'n X ^ 





(aa')^i^^n X -4 In a; 


2 


{aa'^H^ X ^^"^ 


3 


(aaO'i^"^ X ^ 


-9 


{aa'fH^ X -108 Inx 


33 - 18 ln(a) 


{aa')\a'' + a'^)a^/H^ x iff 


(-325+12.^) + 14 i^(^) |ln2(a) 


(aa')'(a2 + a'2)n2 x -iff 


41n(a) + |W(a) 


(aa')3(a2 + a'2)n2x^ 


1 - IHa) 


(aa')='(a2 + a'^)//^^ x 2inx 


(l^)-fln(a) + ilnV) 


+ (>'^)H-U X I 


f - 101n(«) 


(aa')^(a2 + a'2)i/2n x -ISlnx 


-15 +181n(a) 


{aa'f{a^ + a'^)//^ x ^ino; 


(-^) + |Ma)-iln» 


(aa')3(a2 + a'^)//^ x -|| 


-f + f ln(a) 


Total 






Table 5: / (Ext. Operator) x {/(^) - /(V)} each 

contributing nonlocal term. 
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